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We show that photon number measurement can be used to detect superfluidity for a two-band 
Bose-Hubbard model coupled to a cavity field. The atom-photon coupling induces transitions be- 
tween the two internal atomic levels and results in entangled polaritonic states. In the presence 
of a cavity field, we find different photon numbers in the Mott-insulating versus superfluid phases, 
providing a method of distinguishing the atomic phases by photon counting. Furthermore, we ex- 
amine the dynamics of the photon field after a rapid quench to zero atomic hopping by increasing 
the well depth. We find a robust correlation between the field's quench dynamics and the initial 
superfluid order parameter, thereby providing a novel and accurate method of determining the order 
parameter. 



Introduction. - Optical lattices provide very clean and 
highly tunable testing grounds of important strongly cor- 
related Hamiltonians and quantum phase transitions. In 
particular, numerous experimental groups realized the 
bosonic Mott-insulating (MI) to superfluid (SF) phase 
transition in an optical lattice [1] , predicted by the Bose- 
Hubbard Hamiltonian [2J. However, the measurement 
process still stands to be an experimental challenge be- 
cause current methods to probe these phases rely on de- 
structive time-of-flight measurements. To detect super- 
fluidity, experimentalists look for peaks in the resulting 
interference pattern. However, numerical simulations [3] 
have shown that peaks are not conclusive proof of su- 
perfluidty because they may exist even above the critical 
temperature. 

Fortunately, experimentalists can acquire time- 
resolved photon statistics with high precision. Therefore 
one possibility is to circumvent these roadblocks by cou- 
pling the atomic system to light. One could then probe 
the atomic phases by the imprint they leave upon the 
emitted light. Recent advances in engineering strongly 
coupled cavity QED systems provide the opportunity 
to explore fundamental light-matter interactions at the 
quantum level [1H7] . We are encouraged by such exper- 
imental progress to consider harnessing the light-matter 
entanglement to probe novel quantum phases. Here we 
provide a robust nondestructive method to both distin- 
guish the SF and MI phases and determine the order 
parameter. 

To meet our ends, we extend the Bose-Hubbard model 
to include two species of bosons with nearest neighbor 
tunneling. These two species are assumed to correspond 
to internal levels of a single boson, and interaction with 
the photon field induces transitions between the two lev- 
els [Hj. This model can be viewed as an extension of the 
Dicke model, within the rotating wave approximation, to 
bosons in a lattice system [§]. 

As in previous works [T0lfT2] . we find a rich phase di- 
agram for this system, including two component super- 
fluidity and Mott-insulating states. However, previous 



works [ini E] approximate the cavity field to be coher- 
ent, whereas we do not make this assumption. As a re- 
sult, we find differing shapes and positions of the phase 
boundaries when light is present. However, if the field 
is truly a coherent state, the intracavity photon num- 
ber would remain constant after a quench because the 
field and atoms are not entangled. Exact diagonalization 
produces oscillatory behavior, thereby suggesting the co- 
herent field treatment misses a piece of the picture. We 
explore these oscillations by a theoretical analysis that 
allows light-matter entanglement. As we will show, these 
oscillations are a useful probe of the order parameters. 

Our main results are as follows. 1) The average photon 
number is lower in the SF phase than in the MI phase. 
Thus, photon counting is a conclusive method to differ- 
entiate the atomic phase. 2) The entanglement leads to 
nontrivial photon number oscillation after a quench. The 
amplitudes of the oscillations increase with the initial 
value of the atomic tunneling parameter and correspond- 
ingly the initial order parameter. Thus, it is possible 
to deduce the initial order parameter by recording the 
photon number dynamics after the quench. 

Model- Let H a and H b denote the bare Hamiltoni- 
ans of the atomic species. Let Hfi e id denote the en- 
ergy contributions from the photon field and the photon- 
atom coupling. Then, the complete system Hamiltonian, 
assuming negligible cavity dissipation and spontaneous 
emission, is given by: 

H = H a + H b + H field - /iiiVi - fi 2 N 2 (1) 

where, 

H a = -jJ^iata j + h.c.) + e a J2^ (2) 

(ij) i 

H b = -jJ2(4b,+h.c.) + e b J2nt (3) 

(ij) » 

H field = unfil) + -j= 5>la^ + ha\^) (4) 
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FIG. 1. Schematic of Hamiltonian, where the cavity field and 
walls are omitted for clarity. The system consists of bosonic 
atoms with two internal levels. The atoms in their upper level 
are depicted in orange while those in the lower level are in red. 
The cavity field allows transitions between the two internal 
levels, as shown by blue arrows; g denotes the strength of the 
atom-photon coupling. The arrows with parameter strength J 
indicate atomic tunneling to nearest neighbor sites. For sim- 
plicity, we consider only hardcore intra-species interactions 
and zero inter-species interaction. 
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FIG. 2. The phase diagram shows four distinct phases in the 
112/g — J/g plane. Here fi2 controls the number of excitations, 
J is the inter-site tunneling rate, and g is the atom-photon 
coupling. The lower half of the phase diagram has no exci- 
tations, leading to a typical MI-SF transition of uni-species 
hardcore bosons. For larger {12, we predict a field-matter en- 
tangled system. As one tunes J/g in this region, there is a 
MI-SF transition of both atomic species. 



Here, <Zj and bi are the respective ground state and ex- 
cited state bosonic annihilation operators at site i, and tp 
is the annihilation operator of the cavity field. The on- 
site field-atom coupling, -9=(b\a,iijj + bia\%p^), allows for 
atomic transition to the excited state by absorption of a 
photon and the reverse process of emission. The parame- 
ter J is the atomic nearest neighbor tunneling amplitude, 
assumed to be the same for both species. For simplicity, 
we consider hardcore intra-species interactions, and zero 
inter-species interaction. However, our main results do 
not qualitatively depend sensitively upon the strength of 
inter-species interaction. 

Finally, we work in the grand canonical ensemble with 
two chemical potentials corresponding to two conserved 
quantities: the total number of atoms and the total 
number of excitations. The total number of atoms is 
N\ = J^i n i + n i 1 controlled by the chemical potential, 
fii. Following reference [11], the total number of excita- 
tions is N 2 — J2i Vi^i + \{ n \ ~ n t + 1)7 controlled by its 
chemical potential /i2- 

We are primarily interested in the dependence of the 
system on J/g, as we want to focus on the relation be- 
tween the atomic phases and the light-matter coupling. 
Thus, we choose g to be our unit. To keep consis- 
tency and allow comparisons with previous works [X0|. II 1] . 
we also fix the following parameters: fJ-i/g = — .6 and 
5/g = 1, where S = e b — e a — lu. 

To make the atomic tunneling terms more tractable, 
we employ a mean field approximation. Following known 
methods |13l 114] , we arrive at the mean field Hamiltonian 
(Hmf)- 



Hmf = H mf + H b MF + unj>\^ + (5) 

i 

+g(b\a l ip l + bia\i>\) - fixnn - /x 2 n i2 

with 

h mf = e an la - zJ<f> a (a} + a{) + zJ\cj) a \ 2 (6) 
H b MF = e b n ib - zJ4> b (b\ + bi) + zJ\c/) b \ 2 (7) 

where we define a = (a) and 4> b = (b), the respective 
superfiuid order parameters for the lower and upper level 
atoms. 

Note that we also consider xpi to be a local oper- 
ator, rather than global. To justify this, note that 
ip k = J= e lk ri i/>i- However, |k-ri| is negligibly small, 
allowing the replacement ipk — > J^ r . tpi. The following 
sections exhibit the results of our mean field assump- 
tions. However, we note that the trends are in qualitative 
agreement with exact diagonalization results. 

Equilibrium Phases.- We numerically minimize Hmf 
with respect to <j) a and <p b to determine the equilibrium 
phases in the ^2/g — J/g plane. See Figure [2] for the 
phase diagram. In the light-entangled regime, Hmf pre- 
dicts a first order phase transition (Figure [3^,) between a 
MI and a two component SF state. Within mean field, 
the MI state is a lower branch polariton state on each 
site. This is clear from setting <j) a = <p b = in Hmf and 
then noting that Hmf — > J2i-Hi C , where Hf c is sim- 
ply the Jaynes-Cummings Hamiltonian on site i. How- 
ever, the superfiuid phase is a delocalized state of matter, 
thereby necessarily a linear combination of multiple po- 
lariton states. 
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FIG. 3. a) Plot shows how the mean field order parameters of the "a atoms" and "b atoms" vary with zJ/g, indicating a 
first order transition when the system is light entangled, b) Photon number density, (n p h), as a function of hopping strength, 
zJ/g. c) Measure of photon absorption, | (ipi b\ at ) \ / yju^n^na , as a function of hopping strength, zJ/g. Plots a-c indicate that 
the Mott insulating to superfluid phase transition is accompanied by a change in both the intracavity photon number and the 
photon absorption strength. 



Looking at Figure [3}d, we also see that the MI to SF 
phase transition is associated with a discontinuous drop 
in the photon number density, and the number decreases 
continuously as we raise J/g in the SF regime. This be- 
havior results from nonlocal hopping taking precedence 
over the local light-atom interaction with increasing J/g. 
Thus, the photon number tracks the atomic phase, so 
photon number measurement is an easy experimental 
method to accurately detect superfluidity. 

Quench Results- The above features as well as past 
works [151 I16j provide methods to distinguish the two 
atomic phases. However, these works do not distinguish 
the degree of superfluidity. To do so, we need a quan- 
tity that varies sensitively with J. Looking at Figure 
[ik we see that |(^>j&taj}|, a measure of the strength of 
photon absorption, does in fact change sensitively with 
J. Unfortunately, \(ipib\ai)\ is not a directly measurable 
quantity. We therefore consider the equation of motion 
of the intracavity photon number. The number dynamics 
directly depends on ipib\ai according to, 



Our method exploits the quench dynamics to cap- 
ture how the field-matter entanglement affects the pho- 
ton number dynamics. Consider a system in the light- 
entangled superfluid regime. Now consider rapidly rais- 
ing the well depth such that J « 0. The system re- 
tains memory of its initial state, and this state time 
evolves with the zero-hopping Hamiltonian, Hf ina i. We 
phenomenologically include the dissipation through the 
cavity walls, parametrized by k. Thus, the system time 
evolves with the following effective Hamiltonian after the 
quench: 



Hfinai = ^(e a n lQ + e b riib + w^jVi + (H) 

i 



dt 



ig(i>ib\a,i - Wiha\) 



(8) 



Furthermore, the time evolution of ipib\ai depends on 
the initial order parameters, (f> a and (fib- 

d{ ^ ai) = l [H a ^ l b\a^+i[H b ^b\< + - (9) 



d(ipib\ai) 



dt 



= izJ4> a ipib] + ... + -izJ(f>bi>iCLi + ... (10) 



Therefore the absorption strength is the bridge be- 
tween an observable, the photon number dynamics, and 
the desired quantity, the superfluid order parameters. 
Here we propose a way of accessing the information con- 
tained in the photon field's time evolution. 



where k <C g because the system is assumed to be 
strongly coupled. 

Figure [4^, shows the average intracavity photon num- 
ber density as a function of time after the quench. Both 
in the absence of dissipation (red) and in the presence 
of dissipation (blue) we see oscillatory photon number. 
We accounted for dissipation by time evolving with wave 
function Monte Carlo [TTJ [TS] . The probability of detect- 
ing a photon outside the cavity is directly proportional 
to the mean intracavity photon number; thus these in- 
tracavity oscillations are expected to result in oscillatory 
photon counting outside the cavity. 

To understand the oscillatory behavior, we write the 
initial SF state of the system as a linear combination of 
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FIG. 4. a) Quench dynamics of the intracavity photon number density, found from time evolving an initial superfiuid state. The 
red curve assumes zero dissipation. The blue curve accounts for photon leakage by averaging over 10000 wave function Monte 
Carlo realizations (- = .05). The dissipation modifies the photon dynamics by causing an overall decay, but the oscillation is 
preserved for sufficiently small k. b) Fourier transform of the zero-dissipation-dynamics indicates two dominant frequencies. 
The corresponding amplitudes of photon oscillation are denoted A\ and Ai. c) Figure shows the amplitudes, A\ and A2, in 
the presence of dissipation for varying initial hopping. We see that the amplitudes closely correspond to the superfiuid order 
parameters. This trend offers a novel experimental method to deduce the order parameters from photon number dynamics 
alone. Note: for the plots obtained from Monte Carlo, error bars are smaller than the plot line width. 



the eigenstates of Hfi na i: 

|5 imtjai )= Co e^ £ot |0,0,0) + 
+ Cl e^ Elt |0,0,l) + 

—iEkut 



(12) 



£ [ dkt 



fc,i,i) + 



k=l 



+Cfc2e 



+Cfcie lE - * \-,n ph = k) 



-iEf't 



\+,n ph = k) 



where \n p h = k, — ) and \n p h — k, +) are the lower and 

(k) 

upper branch polariton states with k excitations. E_ 

and are their respective energies, and k max is the 
maximum photon number. Furthermore, the Fock states 
( 1 0, 0, 1) and \k, 1, 1)) denote \n p h,nb,n a ). From here, 
(n ph ) (t) is simply: 



(n P h)(t) = A Q 



k=k ma 

E 

fc=i 



A k cos(ft fe t) 



(13) 



with Qk = — E^l' . Therefore each frequency is the 
energy difference between the upper and lower polariton 
states with k excitations. Further, we can see that os- 
cillation occurs only when the initial state has nonzero 
weight in both of these polariton states; more precisely, 
both Cki and Ck2 must be nonzero in order for the number 
to oscillate at frequency f^. For the particular param- 
eters of Figure |4j the initial (mean field) state of the 
system has weight in the following states: \n p h = 1,— ), 



1. 



'ph 



and In 



ph 



There- 



fore, two frequencies are present in the quench dynamics 
(Figure Ep). The dominant frequency corresponds to the 



k = 1 polariton, and the subtler frequency corresponds 
to the k = 2 polariton. 

Next we explore how the dynamics change as the initial 
order parameters change. Consider a system with low 
initial J such that the initial system is Mott insulating. 
As noted in the previous section, the initial state here is 
simply a polariton on each site, \So) — YiiLi \ n ph, — )»• If 
there is zero dissipation, this state is just an eigenstate 
of Hfi na i, so all amplitudes are zero. In the presence of 
dissipation, there is a small offset from zero. 

Now consider starting in the superfiuid regime. At 
smaller initial J, the initial state has more weight concen- 
trated in one particular polariton state over the others. 
That is, from the set of {cki,Ck2} one particular mode 
dominates over the others. For larger initial J, the ini- 
tial state has its weight diffused over multiple polariton 
states, so all the modes from the set {cki, c/^} contribute. 
The validity of these statements is reflected in Figure [3]:, 
where \ {il)ib\a,i) \ decreases with ^ as the nonlocal kinetic 
term increasingly dominates over the local atom-photon 
interaction. Thus, if we start in the superfiuid regime, 
the amplitudes, Af. oc c^ic^, scale with J. 

Figure |4j; exactly exhibits this correspondence between 
the order parameter and the amplitudes of oscillation. 
The amplitudes, A\ and A 2l are extracted by fitting the 
Monte Carlo curves to the following: 

(np h ) (i) = A - (A 1 cos(O x f) + A 2 cos(fV))e" c * (14) 

We emphasize that Figure|3J: relates to the dynamics of 
the intracavity photon number density. The total number 
within the cavity scales with the number of lattice sites. 
As a result, the amplitude trends in Figure [4]: would be 
far more pronounced in actual experiment. 
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Thus, we come to the central conclusion of this section; 
the correspondence between the initial order parameters 
and the photon number quench dynamics allows elucida- 
tion of the initial quantum phase from photon statistics 
alone. 

Experimental Method to Determine Superfluid Order 
Parameters. - Our theoretical analysis indicates that pho- 
ton number measurement allows detection of superflu- 
idity in a two band Bose Hubbard system coupled to 
light. Furthermore, the light-matter interaction bridges 
the photon number dynamics and the superfluid order 
parameters. As a result, measuring the photon number 
amplitudes of oscillation after a quench provides an easy 
experimental measure of the order parameter. Here we 
note that we expect realization of our calculations to be 
possible within the near future. In particular, the group 
of Thompson and collaborators [7] succeeded in placing 
an optical lattice in a cavity, albeit a thermal system. 

Future directions for this work can treat dissipative 
effects more carefully. For our calculations we assumed 
~ = 20 and negligible decay of the upper state, 7. Ex- 
periments may benefit from knowing how 7 and k affect 
the above trends. 
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